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Abstract 

In this work, we study the Benjamin-Bona-Mahony hke equations with a fully nonlinear dispersive term 
by means of the factorization technique. In this way we find the travelling wave solutions of this equation 
in terms of the Weierstrass function and its degenerated trigonometric and hyperbolic forms. Then, we ob- 
tain the pattern of periodic, solitary, compacton and kink-like solutions. We give also the Lagrangian and 
the Hamiltonian, which are linked to the factorization, for the nonlinear second order ordinary differential 
equations associated to the travelling wave equations. 

Email:kuru @ science.ankara.edu.tr 



1 Introduction 

The nonlinear dispersive equations K{m, n) (or generalized Korteweg-de Vires equation) 

ut + iu"')x + {u"-)^^^ = 0, m > 0, 1< n < 3 (1) 

have been dealt with by Rosenau P. and Hyman J.M. (|T|. They built travelling wave solutions of ([T]i with 
compact support and called them compactons. Such compactons are solitary waves but, besides the compact 
support, they have some different features from the solitons. So, while the width of the soliton depends on its 
amplitude, the width and amplitude of a compacton are independent. However, compactons sometimes behave 
as solitons: they also collide elastically and their amplitudes depend on speed of the wave. In the literature 
there are many studies on K{m, n) equations and their compacton solutions H] |2] [3] |4] |5] IS [2l Ej . 

The Benjamin-Bona-Mahony (BBM) equation has a higher order nonlinearity of the form 

Ut + Ux + a Ux - u^.j,t = 0, n > 1 (2) 

where a is constant |'9',T^. This equation is an alternative to the Korteweg-de Vries equation, and describes the 
unidirectional propagation of small-amplitude long waves on the surface of water in a channel, hydromagnetics 
waves in cold plasma and acoustic waves in anharmonic crystals. Thus, in this sense, it has some advantages 
compared with the Korteweg-de Vries equation. 

Now, we will consider the BBM-like equations, B{m^ n), with a fully nonlinear dispersive term, similar to 
the K{m, n) equations, 

wt + a(u™)^ - (u"):rxt = 0, TO, n>l. (3) 

There are some studies related with these BBM-like equations lfT3lfT4l[T5l[T6ll . in general using the sine-cosine 
and the tanh method, in order to get the travelling wave solutions of this family. In this respect, let us remark that 
in the literature there are many methods to obtain exact or special travelling wave solutions (soliton solutions) 
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of nonlinear equations ifTTl . but there is not a standard method. So, it is important to investigate new methods to 
get solutions in different ways. In this work, we deal with BBM-like equations systematically by means of the 
factorization technique ifTOl [iTl [TSl [T9I to find the travelling wave solutions. When looking for the travelling 
wave solutions, the BBM-like equations will reduce to second order nonlinear ordinary differential equations 
(ODE). Then, we can apply straightforwardly the factorization to these ODE's. 

In section 2, we briefly introduce the factorization technique for a class of nonlinear second order ODE's. 
Then, in section 3 we apply this method to the BBM-like equations. In section 4, we obtain the travelling 
wave solutions of this equation in terms of Weierstrass functions which cover the pattern of solitary wave 
and periodic solutions. We also construct compactons and kink-like solutions for this family and we give the 
graphics of these solutions in Fig. [T] (periodic). Fig. |2] (compactons). Fig. [3] (kink-like) and Fig. |4] (solitary 
wave pattern). We give the Lagrangian and the Hamiltonian corresponding to the second order nonlinear ODE 
and their relation with the factorization in section 5. Finally, in the last section we end this work with some 
conclusions. 



2 Factorization of the nonlinear second order ODE 



Let us consider, the nonlinear second order ODE of the special form 

(PW dW , , 



(4) 



where /3 is constant and F^W) is an arbitrary function of W . The factorized form of this equation can be 
written as 



dQ 



dO 



h{w,e) 



w{e) = 0. 



(5) 



Here, /i and /2 are unknown functions that may depend explicitly on W and 6. Expanding (O and comparing 
with dill, we obtain the following consistency conditions 



F{W) dh . ^ d{Wh) _ 
If we solve (|6|l for /i or /2, it wiU aUow us to write a compatible first order ODE 



W(9) = 



(6) 



(7) 



that provides a solution for the nonlinear ODE (|4]l ifTOl [iTl [TSl [191 . In the applications of this paper /i and /2 
will depend only on W. 



3 Factorization of the second order nonUnear ODE's corresponding to 
the BBM-like equations 

Let us assume that ^ has travelling wave solutions 

u{x,t) ^ (f){^), £,=hx + wt (8) 

where h and w are real constants. Substituting ^ into ^ with the condition m = n and after integrating, we 
obtain the second order nonlinear ODE 

(<^")45- A0-B0" + i? = O (9) 
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together with the constants 

n.^ w nw h'^ w 

here R is integration constant. Now, if we introduce the following natural transformation of the dependent 
variable 

= w{e), i = e (11) 

then, Eq. ^ becomes 

'i^-AW^''''~BW + D = Q. (12) 
Next, we will apply the factorization technique to Eq. (fT2] l. In this case, (3 — Q and F{W) is 

i^(M^) = -(AM^i/" + BVK-i:>) (13) 
therefore, the consistency conditions take the form: 

/1/2 = -A Vl/(i"")/" - B + D (14) 

= (15) 
Substituting ( fT4] i into (flST l. we get only one first order consistency equation 



whose solution gives the unknown function /i. 



fl + h - A 1^(1-")/" -B + D W-^ = (16) 

dW 



fi(W) ^±^JbW^ -2DW +^^W'-''+^y'' + C (17) 

^ ^ V n + 1 

where C is an integration constant. Then, we can write the first order ODE (|7|i 



\ BW^-2DW+^^W^"+^)/'' + C = 0. (18) 
de \ n+l 

In order to solve this equation for 14^ in a more general way, let us take W in the form W — cp^ , p 0,1 so, 
the first order ODE ( fTsT i can be rewritten in terms of ip as 

.dip.y B o 2D o_„ 2nA C n_o„ 

aO (n + ljp p 

If we want to guarantee the integrability of ( fT9] l. the powers of have to be integer numbers between and 4 
II20I . Having in mind the conditions on n {n > 1) and p{p ^ 0, 1), we have the following possible cases: 

• If C = 0, D — 0, then we can choose p in the following way 

. 2n n n 2n . 

,-1 ,1 ,1 }■ (20) 

1 — n 1 — n 1 — n 1 — n 

Here, we deal with only the case p — —j^, since the cases p — and p = give rise to the 

same solutions for (|3]l. Thus, Eq. ( fT9] l takes the form forp = —j^' 

ctb' 71^ n (n + 1) 

Notice that for the special values of n = 2, 3, we have the B{2, 2) and _B(3, 3) equations and the obtained 
solutions can be read as particular cases of this situation. 

• If C = and D 7^ 0, we have the special cases: [p — 2, n — 2) and [p — —2, n ~ 2) which correspond 
to the 5(2,2) equations. Since these two cases also give rise to the same solutions for we will 
consider only one of them. Now, Eq. ( fT9l ) has the following form for [p — 2, n — 2) 

.dip. 2 B 2 A D 
(J) =4^ +3^- 2- 
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4 Travelling wave solutions for BBM-like equations 



In this section we will obtain the solutions of the differential equations (ISTT i and ( |22] | which allow us to get the 
travelling wave solutions of B{n, n) (O. To solve this type of differential equations, it will be useful to recall 
some properties of the Weierstrass functions 11211 l22ll . 

Let us consider a differential equation with a quartic polynomial 

{-^Y = P{^) = flo'^'' + 4ai (^^ + 6a2 </3^ + 4a3 V9 + 04 . (23) 
^ ad ' 

The solution of this equation can be written in terms of the Weierstrass function p{0] 32 , 53) where the invariants 
g2 and 33 of ( |23] ) are 

32 = flo 04 ~ 4 ai 03 + 3 53 = oq 02 04 + 2 ai 02 03 — — ao 03 — a\ 04. (24) 
Then, the solution ip can be found as 

V{0) = V'o + Ji^'lv'o) (p(e;52,53) - ^^"(^0)) (25) 

where the prime (') denotes the derivative with respect to p, and ifQ is one of the roots of the polynomial 
P{f) ( 1231 ). The discriminant A = ~ 27 = 0, allows us to express Weierstrass functions in terms of 
trigonometric and hyperbolic functions 11211 l22l 

p(6l; 12 6^ -8 6^) = 6 + 3 bsinh^^p ^)i/2^] (26) 

p(6'; 12 8 6^) = -6 + 3 6sin-2[(3 6)i/2g)], (27) 

Now, we will examine by separate each of the two cases considered in the above section. 

4.1 C = 0, D = 

In this case we have second order polynomial 



n (n + 1) 



p{^) = '" ;^ ^ ^^+ ^^' (28) 



whose roots are 

(^0 = 0, ^0 = -^^^^. (29) 
[n + l)i( 

The invariants have the form 

^ 12 n4 ' ^ ^I^;^ ^3^) 

and discriminant A 0. 

Then, we can find the solutions from (|25] | in terms of the Weierstrass functions: for (po = 

^n(n-\f A 

'^^ ^ - " (n + 1) \[n - 1)2 B - 12 n2 p[Q- 52, 53)] ^ ^ 

and for (^0 = -(^^f^s 

"2(n-l)2S + 12n2p(0;52,53) 



2 71 A 



(n-l)2B-12n2p(0;g2,g3) 



(32) 
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Since the discriminant A = 0, we can substitute (|26] | and ( l27l i into the solutions (l3Tl i and ( l32b . leading to 
the following type solutions of (l2Tl i: a) periodic solutions 



Vi0) = 

for B < 0, and b) hyperbolic solutions 



sm 



2nA 
{n + l)B 

2nA 
(n+ 1)5 



2nA . , 2 
Sinn 



1) 



(n + 1) B 
2nA 



cosh' 



2n 






-1) 


\ 2n 




fy/B{n- 




{ 2n 









2n 



for i? > 0. Notice that the first two real solutions (|33] | and ( |34| | are simply related by a transformation. 
Now, taking into account (O, ( fTTT i and W — ip^ , p — n/{n ~ 1), the solution of (|3]l can be written as 

Substituting ( l33T l and ( |34] | into ( |37] |. we obtain periodic solutions of the equation (O 

l/n-l 



(33) 
(34) 

(35) 
(36) 



u{x, t) 
u{x, t) 



2nA 
' {n+l)B ' 

2nA 

{n + l)B 



cos 



{n - 1) 
2n 

^{n- 1) 
2n 



(/i a; + w i) 



{hx + wt) 



l/n-l 



(37) 



(38) 



(39) 



for B < 0. Then, combining the trivial solution u{x, t) = Q with (l38b and ( [39l ). we have the compact support 
solutions for ([3]l in the following way IT], 



u{x, t) = 



u{x, t) = 



-^sin^{^^^^ihx + wt)) 



2nA 

{n+l)B 



{hx -\- wt) 



1/ri-l 



l/n-l 



0< + < 

otherwise 

otherwise. 



(40) 



(41) 



For the trivial solution u{x, t) = (— we also get the compactons: 



u{x, t) 



u{x, t) = 



( 2nA \l/n-l 
y (n+DB' ' 



l/n-l 



2rLi4_ ^Qg2 



(n+l)S 

2n/ 
(n+1) 



l/n-l 



/ 2nA \l/n-l 
y (n+DB' ' 



\^{hx + wt)\<j^ 
otherwise 

0<^{hx + wt) 
otherwise. 



(42) 



(43) 
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Now, we can construct kink-like solutions ll23l . combining the non trivial solutions dSSl l and ( |39] ) with both 
(different) constant solutions, u{x, t) = and u{x, t) = (— )^^"^"^. as follows 



r 0, 



n+l)_B \^ 



(n+ 

2n 
' (ri+1) S 



(k-1) 
2n 



(/i a; + w i) 



l/ri-l 



/ 2nA M/ra-1 



t) — < 



2„A ,^Qg2 



(n+1) B 



l/n-l 



0, 



{hx + wt) < 

< V^ihx + wt) < 
^{hx + wt) > ^ 
^/^{hx + wt) < 

0<V^(hx + wt) < 
V^{hx + wt) >J^y 



(44) 



(45) 



When we put i35[ and ( l36l l in dJTl l. we have the solitary wave pattern solutions of hyperbolic type for 
In order to get real solutions, it is necessary to examine the solutions for n even and n odd. Then, for n even 
and B > we have the solutions 



u{x, t) 



u{x, t) 



2nA . , 2 / \/^("- 1)., 

■sinh^ I i -{hx + wt) 



{n + l)B 



2n 



2nA ,2//B("-1)/, 

■cosh^ I ^ -{hx + wt) 



in+l)B 



2n 



l/n-l 



1/ri-l 



(46) 



(47) 



However, if n is odd, we have the solution ( l46b provided yl > 0, while the solution (07]) is valid only when 
^ < 0. 

The solutions for the B{2, 2) and the -6(3, 3) can be easily read from the above solutions of (O. 

u{e) 




-20 -10 



10 20 




-20 -10 



Figure 1: Plot of trigonometric solutions tSm and (l39t forh= 1, w — 1, a = —1, n — 2 



4.2 C = 0, D^O 

The case n = 2 corresponds to the -8(2, 2) equation and, according to ( l22b we have another kind of solutions 
for this equation. In this case, the polynomial 

P{^)-^v' + ^^-§ (48) 

has two roots 



± -A±ViA^ + 18BD 
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Figure 2: Plot of compactons ( l40l l and (l42l l for /i=l,?i; = l,a = — 1, ?i = 2. The plot of compactons (HTt 
and (|43] | are translations of ( l40l i and (|42] |. respectively. 




u(e) 




-20 -10 



10 20 



-20 -10 



10 20 



Figure 3: Plot of kink-like solutions ( |44] | and i45[ for /i=l,w = l,a = — l,n = 2. 



From (l24l l. the invariants for Weierstrass functions take the values 



92 



93 



192' 13824 
and the discriminant also vanishes A = 0. Then, from dZSl ) the solutions read 

2A VaWTTsWd; f 5 B + A8 p{0; g2, 93) 



<p± = T 



3B 



B - 48 p(0; 52,53) 



(50) 



(51) 



Since the discriminant equals zero, we can express these solutions in terms of trigonometric and hyperbolic 
functions using the relations (|26] | and (|27] |: 



for _B < and 



2 A V4A2 + 18BD 
"35* 35 — 



2A V4A2 + 1851) ,,V5 
'35± 35 — 



(52) 



(53) 



for 5 > 0. Hence, from (|37] | we have ^(x, i) = (/>(^) = (^(0), 9 = ^ = hx + wt, and the solitary wave pattern 
and trigonometric solutions of (O can be found substituting ( |52] | and ( |53] | into this expression. Here, again 
combining the trivial solution u{x,t) = — with the trigonometric solution, we can construct compactons as 



u^{x, t) = 



± ^'a^i^ — cos(^ {hx + wt)) , \V^ihx + wt)\ <TT 



2 A 
3B ' 



(54) 



otherwise . 
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u(e) 



-3 -2 -1 




1 2 3 

Figure 4: Plot of solitary wave pattern solutions (l46T l and (|47] | for /i = 1, w = 1, a = 1, n = 2. 



Kink-like solutions can also be obtained by taking into account another trivial solution u{x,t) 
^i^'+iSBD together with u{x,t) ^ - 



2 A . 
3B- 



2 A 

' 3B' 



'B (hx + wt) < -TT 



u'^ix, t) = < 



_ 2_A _|_ ^iA^ + lSB D 
3B 3B ' 



(/l X + W t) > 

where AA^-\-l%BD > and -B < 0. The solitary wave (hyperbolic type) pattern solutions are 



3B 



cosh(— ^ (ft-a; + 



with the conditions 4 + 18 B D > and B > 0. 



2 A 

"3S 



± 



_2|^ Vimm cos(^(/ix + wt))l , < \/^(/i.T + w<) < (55) 



(56) 



5 Lagrangian and Hamiltonian 

From the second order nonlinear ODE ( ITS] ), the Lagrangian can be written as 

Lw{W, VKe, 6*) = i (^"^ + 1^"^ + B _ 2 £) ) . (57) 

Then, the canonical momentum is Pw ~ ^ewl ~ Hamiltonian Hw = WePw — Lw has the form 



1 / n 2nA n + l „ 

"2 y^r— - BW^ + 2DW 



Hw{W,Pw,0)^t;[Pw 

2 \ n + l 

Since Hw does not depend on the variable 0, it is a constant of motion Hw = E 

2 \ \de J n+l ) 



(58) 



(59) 



This equation also gives rise to the first order ODE ( fTSl l with the identification C = 2 E. We can express this 
constant of motion Hw — E as a product of two independent constants of motion 



where 



2nA 
n+l 



+ BW^ ~2DW) e^(^) 



(60) 
(61) 
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h = [We + J ^ +BW^~2DW) e'^^'^ 




(62) 



and S{d) has the form such that Ii and I2 satisfy dIj/dO = 0, j = 1, 2 




AW- +BW 




+ BW^ - 2DW 



6 Conclusions 

In this work, we have investigated the travelUng wave solutions of the BMM-Hke equations by means of the 
factorization technique. We have factorized the nonhnear second order ODE's corresponding to the BMM-Uke 
equations. Then, we have obtained the solutions in terms of Weierstrass functions giving rise trigonometric 
(periodic) and hyperbolic type solutions (solitary wave pattern) and we have constructed compactons and kink- 
Uke solutions for the BBM-like equations. In addition to these we give the Lagrangian and the Hamiltonian 
corresponding to the second order ODE. We have seen that the first order ODE which gives rise to the solutions, 
can also be obtained from the constant of motion E corresponding to the Hamiltonian. Then, we have written 
this constant of motion E as a product of two independent constants of motion. We note that this technique is 
more systematic than others previously used for the analysis of this type equations. In general, we have more 
general solutions and we have recovered all the solutions previously reported lfT4l [TSl [161 . Finally, we must 
mention that this method can be easily implemented to the other nonlinear equations, in particular we plan to 
give some results on the B{m, n) equations with m 7^ n in a future publication. 
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